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IMPLICITIZATION OF RATIONAL HYPERSURFACES VIA 
LINEAR SYZYGIES: A PRACTICAL OVERVIEW 

NICOLAS BOTBOL AND ALICIA DICKENSTEIN 


Abstract. We unveil in concrete terms the general machinery of the syzygy- 
based algorithms for the implicitization of rational surfaces in terms of the 
monomials in the polynomials defining the parametrization, following and ex¬ 
panding our joint article with M. Dohm. These algebraic techniques, based on 
the theory of approximation complexes due to J. Herzog, A, Simis and W. Vas- 
concelos, were introduced for the implicitization problem by J.-P. Jouanolou, 
L. Buse, and M. Chardin. Their work was inspired by the practical method of 
moving curves, proposed by T. Sederberg and F. Chen, translated into the lan¬ 
guage of syzygies by D. Cox. Our aim is to express the theoretical results and 
resulting algorithms into very concrete terms, avoiding the use of the advanced 
homological commutative algebra tools which are needed for their proofs. 


1. Introduction 


Let K be a field. We can assume K = Q (or any computable field) when dealing 
with implementations. All the varieties, rings and vector spaces we will consider 
are understood to be taken over K. Consider a rational parametrization 


( 1 . 1 ) 


/ 


S=(S1,S2) 


( fijs) / 2 (s) /a(a) A 

\fo{sy MsY Ms)) 


of a (hyper)surface y := {F = 0) C A^, where F G K[ri,r 2 ,T 3 ] is a non-constant 
polynomial and fi G IK[si, S 2 ]. (As usual, the dashed arrow means that / is defined 
on a dense open set of K^.) An important problem in computer aided geometric 
design is to switch from parametric to implicit representations of rational surfaces 
[Hof89] . that is the parametrization f is assumed to be known and one seeks for 
the implicit equation F (which is defined only up to multiplicative constant). In 
fact, we will assume that / is given and our aim will not be to get the implicit 
equation F oi written in terms of its monomials, but a matrix representation of 
the surface. 


Definition 1.1. A matrix representation M of is a matrix with entries in 
K[ri,T 2 ,T 3 ], generically of full rank, which verifies the following condition: for 
any point p G the rank of M{p) drops if and only if p lies on SA. 

The use of matrix representations goes back to Manocha and Canny |MC91j . 
and to Chionh and Goldman |CC92) . Having the matrix M is sufficiently good for 
many purposes and it is cheaper to compute. The well-developed theory and tools 
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of linear algebra can be applied to solve geometric problems. We can certainly use 
the (numerical) rank dropping condition in Definition o to check membership in 

and, moreover, the whole structure of minors of M is related to the singularities 
of the parametrization |BBC14j and gives a way to invert it when the fiber has a 
single point |Busl41 IBBC14| . Matrix representations are also useful for solving 
intersection problems as is shown in |ACGVS(T7l ITBM091 IDTFGVS13I IBusl4| . 
Much of the computational difficulty in these problems lies on computing ranks 
for polynomial matrices (cf. |HS99| as well as Section 5 in the nice and interesting 
paper IBusldj l. 

The motivation for this paper is to present in the simplest possible terms proce¬ 
dures for the implicitization of rational surfaces via matrix representations, based on 
the syzygies ..., / 13 ) of the input polynomials, that is, 4-tuples of polynomials 
in the s variables verifying the linear relation X]i=o theoretical jus¬ 

tification is not naive and requires a good command of techniques of (homological) 
commutative algebra. However, the algorithms do not require a heavy background 
and are easy to explain. We will show that they perform very well, and moreover, 
they work even better in the presence of base points. 

Call Ti,T 2 ,T 3 the coordinates in the target of /. Our question is an instance 
of elimination of variables, where we want to find the algebraic relations among 
the variables Ti,T 2 ,T 3 under the assumption that /o(s)7i — fi{s) = 0,z = 1,2,3, 
for some s in the domain of /. The eliminant polynomial by excellence is the 
determinant det(A), a polynomial with integer coordinates on the coefficients of 
a square matrix A, which vanishes on those coefficients for which there exists a 
nonzero solution x to the equations A ■ x = 0. Elimination of variables is done in 
the literature through different incarnations of the following general strategy: 

(1) Reduce the problem to a linear algebra problem. 

(2) Hide the variables one wants to eliminate in the (typically monomial) bases. 

(3) Use determinants. 

This strategy is also the core in our syzygy-based algorithms. 

The following short account of the approach of the use of syzygies in our context 
is reconstructed from David Cox’s lecture at the Conference PASI on Commuta¬ 
tive Algebra and its connections to Geometry honoring Wolmer Vasconcelos, held 
in Brazil in 2009 |CP111 Mini Course 1]. The use of syzygies for the impliciti¬ 
zation of (conic) surfaces goes back to Steiner in 1832 |Ste32| . In 1887, Meyer 
describes in |Mey87| syzygies of three polynomials and makes a general conjecture 
proved by Hilbert in 1890 |Hil90| . Surface implicitization by eliminating param¬ 
eters was studied by Salmon in 1862 |Sal58) and Dixon in 1908 using resultants 
[Dix08) . In 1995, Sederberg and Chen reintroduced the use of syzygies, by a 
method termed as Moving curves and surfaces |SC95| . Cox realized they were 
using syzygies |Cox01j . and produced several papers with other coauthors (Buse, 
Chen, D’Andrea, Goldman, Sederberg, Zhang [BCD03[|Cox03a[ICSC98llZSCC03) b 
In 2002, Jouanolou and Buse |BJ03| abstracted and generalized on a sound basis 
the method of Sederberg-Chen via approximation complexes, a tool in homological 
commutative algebra that had been developed by Herzog, Simis and Vasconcelos 
[HSV83bllHSV82IIHSV83aj . Buse, Chardin, Jouanolou and Simis produced further 
advances in the homogeneous case |BjfmiBho^l(ffia06llBrnomiBCS10| . Goldman 
et al. studied the cases of planar and space curves |JG09[ IHWJGIOI IJWG10| . A 
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generalization of the linear syzygy method when the support of the input polynomi¬ 
als is a square (that is, bihomogeneous of degree (d, d)) was proposed by Buse and 
Dohm |BD07| . and for any polygon by Botbol, Dickenstein and Dohm |BDD09j . 
and Botbol |Botn9LlBHmillBotllb| . This method is particularly adapted when 
the polynomials defining the parametrization are sparse, which is often the case. 
This will be our point of view in this article. So, we want to solve the following 
problem. 


Problem. Given a rational parametrization f as m(l.l|, find a matrix represen¬ 
tation M of the surface 5^ by means of syzygies and the monomial structure of 
/o, • • •,/a- 


The main general algorithmic answer to this problem is given in Algorithm |3.1| 
(see Theorem 3.3). Our assumption that the dimension of is 2 is equivalent to 
the fact, when we extend the map to the algebraic closure K of K, that for almost 
all p = f{s) in the image of /, the number of preimages by / is finite. This number 
is called the degree of / and noted deg(/). The matrix representations M oi 
provided by Algorithms |3.1| and |3.6| moreover satisfy that the greatest common 
divisor of all minors of M of maximal size equals ^ 


We present in Section the first naive linear algebra algorithm to compute the 
implicit equation A, which requires to solve a huge linear algebra system. Moreover, 
this naive method “forgets” the parametrization and thus in general it is not useful 
in Computer Aided Geometric Design. In Section we recall previous results on 
the implicitization of curves and surfaces using syzygies and present our general 
methods of implicitization via linear syzygies, which requires to solve considerably 
smaller systems. We highlight in Section [3.1| the main elimination step, which was 
termed instant elimination in [Eis04) (see also the references therein). 

In Section |3.3| we present in Theorem |3.8| a rehnement of Theorem |3.3| for bi¬ 
homogeneous parametrizations, in the same spirit. Technicalities are avoided in 
our presentation in these sections, and in particular in the statement of our main 
results Theorems 13.31 and 13.81 

Detailed hypotheses and proofs are deferred to Section where we introduce 
the necessary background on toric geometry. We collect in Appendix a general 
overview of the rationale of the tools and results from homological commutative 
algebra required for the proofs. A reader only interested in the application of our 
results, can skip these two sections. 

Section illustrates the practicality and advantages of our approach. For our 
computations, we use implementations in Macaulay 2, which need different type of 
homogenizations to use current routines (via a toric embedding or a multihomoge¬ 
nization via an abstract toric Cox ring) [BDlOl IBotlO| ^ For the best performance 
of our algorithms, it would be important to design ad-hoc structured linear algebra 
strategies to compute syzygies in the sparse case. 


2. A NAIVE LINEAR ALGEBRA ANSWER 

The convex hull in K" of the exponents of the monomials occurring in a non zero 
(Laurent) polynomial h in n variables is called the Newton polytope Af{h) of h. 
When h is a polynomial in (si, S 2 ) of degree (at most) d, its Newton polygon A/’(F) 


^Routine updates at: http://mat6.din.uba.ar/~nbotbol/Macaulay2/BigradedImplicit.m2, 
http: //mate. dm.uba. ar/~nbotbol/Macaulay2/MatrixRepToric .m2. 
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is contained in the triangle Ad with vertices (0,0), (d, 0), (0, d). The Euclidean area 
vol(A(i) of this triangle is d‘^/2 and its lattice area volz(Ad) is equal to 2-d‘^ j2 = d^, 
which is always an integer. 

We have the following classical result (c.f. for instance |BJ03p : 

Theorem 2.1. For generic polynomials /o ,...,/3 of degree d, the degree of the 
implicit equation F is df and its Newton polytope is the tetrahedron with vertices 
(0,0,0),(d2,0,0),(0,d2,0),(0,0,d2). 

In the sparse case, the following generalization holds |SY94] . 

Theorem 2.2. For generic polynomials /g, ■ ■ ■, /a with the same Newton polygon 
P, the degree of F is the lattice area v = volz(P) and its Newton polytope is the 
tetrahedron with vertices (0,0, 0), (i;, 0,0), (0, u, 0), (0,0,1;). 

A first naive algorithm would then be the following. Assume the Newton poly¬ 
tope Af (F) of F is known (as in the previous theorems) and number mi,..., mA? € 
the integer points (also called lattice points) in Af{F). Consider indeterminates 
c = (ci,... ,cjv) and write F = CiT^' . Substitute T = f{s) and equate to 
zero the coefficient of each power of (si, S 2 ) that occurs (clearing the denominator). 
This sets a system C of linear equations in c, with solution space of dimension 1. 
Any nonzero solution c of T will give a choice of implicit equation F. 

This solves the problem, but, which is the size of this linear system £? 

The number of lattice points in A^ equals . In the the sparse case, the 

number of lattice points of a given lattice polygon P can be computed via a theorem 
of Ehrhart valid for any n [Ehr67| . which amounts to Pick’s formula in the case 
n = 2. Given a positive integer t, we denote by tP the Minkowski sum of P with 
itself t times, i.e. tP — {pi + ■ ■ ■ + pt, Pi & P for i = 1,... ,t}- The number of 
lattice points in tP equals 

(2.1) #(tP n Z2) = Yo\iP)f + ^ voh{dP)t + 1, 

where volz(i9P) denotes the number of lattice points in the boundary of P. In 
particular, #(P n Z^) = vol(P) -b ^ volz(9P) -b 1. 

The proof of the following result is straightforward: 

Lemma 2.3. In case /o ,...,/3 are generic polynomials of degree d in (si,S2), 
the number of unknowns in the linear system C in the coefficients of the implicit 
equation F is ('^ (approximately d^/6) and the number of equations is (f^ 2 *"^) 
(approximately df 12). 

For any lattice polygon P and generic polynomials fi with Newton polytope P, the 
linear system C has (approximately volziP)^/Q) variables and bbMPL_|_ 

voiz^p) _|_ I equations (approximately yoIz{P)^(2). 

We will see in Remark [3.5| of Section|^that the size of the involved linear systems 
in the syzygy based methods is drastically smaller. 

3. The main algorithm based on linear syzygies 

Our main result is Theorem |3.3[ which has a wide applicability. We distill and 
state it in naive terms, which do not call upon the more sophisticate tools recalled 
in Section and Appendix required for its proof. This is why we postpone the 
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detail of Hypotheses |4.4| and |4.8| until Section]^ Our approach is an inhomogeneous 
translation of the basic general algorithm for the sparse case in |BDD09| , which were 
inspired by the methods |BJ03j for classical homogeneous polynomials. 

Before moving to the implicitization of rational surfaces, we recall the practical 
approach of moving lines proposed by Sederberg and Chen |SC95| for the implicit¬ 
ization of planar curves. 


3.1. Curves. A planar rational curve ^ over a field K is given as the image of a 
map 


f fljs) f2{s) \ 

KfoisY Ms)]’ 

with fi G K[s] polynomials of degree d in s. We can assume without loss of gen¬ 
erality that gcd(/o, /i, / 2 ) = 1. Remark that a linear syzygy can be represented as 
a linear form L = h^To -|- hiTi + h 2 T 2 in the new variables T = {Tq,Ti,T 2 ) with 
hi G ]K[s] such that 

^ h^f^ = 0. 

i=0,l,2 

With this incarnation, a linear syzygy was termed a moving line. For any G N, 
consider the finite-dimensional K-vector space Syz{f)^ of linear syzygies satisfying 
deg(hj) < iz, and call N{iz) its dimension. 



Pick a K-basis h® = (/iq, ..., h\),i = 1,..., N{iz) of Syz(/)i^. Consider the mono¬ 
mial basis {1, s,..., of polynomials in s of degree at most v and write for each 
syzygy hh 


L. = L.(s,T)= ^ h]is)T,= Y. 

i-0,1,2 \k^0 / 


"EE 4i.t, I 


fc=0 \i=0.1,2 


Let Mu be the N{i/) x (zz -|- 1) matrix of coefficients of the Li’s with respect to the 
basis {1, s,..., s''}: 


M, 





Observe that the variable s has disappeared. This is the main elimination step! 


It is known that for iz > d— 1, the matrix M^, is a square matrix with det(Mi^) = 
7 cdeg(/)^ where F is an implicit equation of In case v > d, then Mi, is a non¬ 
square matrix with more columns than rows, but still the greatest common divisor 
of its minors of maximal size equals both cases, for zz > d — 1, a point 

P G lies on iff the rank of M„(P) drops. 

In other words, one can always represent the curve as a square matrix of linear 
syzygies, which gives a matrix representation of the implicit equation. In principle, 
one could now actually calculate the implicit equation, but the matrix M^, is easier 
to get and well suited for numerical methods |ACGVS07] . As we remarked in the 
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surface case, testing whether a point p lies on the curve only requires computing 
the rank of M^, evaluated in p. Also, the singularities of C can be read off from M^, 

In the absence of common zeros of /o,/i,/ 2 , it is possible to find the implicit 
equation via a resultant computation. Note that for a parametrization with poly¬ 
nomials of degree d, the Sylvester resultant matrix uses a matrix of size 2d, while 
the syzygy method uses 2 matrices of size d, as the Bezout resultant. 


3.2. The general method of implicitization via linear syzygies for surfaces. 

Assume we are given a rational parametrization of a surface .5^ as in (1.1). We aim 
at finding a matrix representation for Note that we can in principle assume 
that (/o,..., /a) are Laurent polynomials admitting negative exponents, but after 
multiplying them by a common monomial, we get a new rational parametrization 
of =5^ defined by polynomials fi G K[si,S 2 ]- We will then assume, without loss of 
generality, that we are have a lattice polygon P which lies in the first orthant of 
and contains the Newton polytopes of /o, ■ ■ ■, /a- 

We saw that in the curve case, it is always possible to find a square matrix 
representation. In the surface case, however, linear syzygies provide in general 
rectangular matrix representations and the implicit equation (raised to the degree 
of the map /) equals the great common divisor of the maximal minors (or the deter¬ 
minant of a complex). A recent paper by Buse |Busl4] presents a very interesting 
square matrix representation out of a matrix representation M when we work over 
the real numbers, by considering the square matrix MM*. This approach is natural 
because of the properties of the rank of a real matrix with respect to its singular 
value decomposition. The determinant of MM* gives an implicit equation for 
(in general, it gives F with multiplicity), which is moreover a sum of squares. As 
Buse observes, for complex matrices it would be enough to replace the transpose 
M* by the conjugate transpose. 

The use of quadratic relations (i.e. linear syzygies among the products fifj of 
any of two of the polynomials fi defining the parametrization) was proposed to 
construct square matrices |SC951 ICSC981 ICoxOll ID’AOll lAHWOS] . Khetan and 
D’Andrea generalized in 2006 [KD06j the method of moving quadrics to the toric 
case. The choice of the quadratic syzygies is in general not canonical and the cost of 
computing syzygies is increased. Note that syzygies in (/q, ..., /s) and the implicit 
equation F have a common shape. Indeed, linear syzygies h = (ho ,..., / 13 ) of degree 
1 / correspond to polynomials H(s, T) = J2i=o such that = Oj 

with deg(iL) in the s variables equal to iz, and deg(iL) in the T variables equal to 
1. Also, quadratic syzygies of degree v' correspond to polynomials H(s,T) = 
J2t<j=ohijis)TiTj such that deg(iJ) in the s vari¬ 

ables equal to F, and deg(iL) in the T variables equal to 2. The implicit equation 
(of degree D) is a polynomial F[{s, T) = X]|ct|<£> such that h^fi* (s) = 0. 
Thus, deg(iL) in the s variables equals 0, and deg(iL) in T variables equals D. So 
to go from linear syzygies to the implicit equation, in some sense one has to play 
the game of lowering the degree in the s variables to 0 (which increases the degree 
in the T variables up to D). 

We now present our main general algorithm to construct matrix representations 
of parametrized surfaces. Clearly, given any lattice polygon P C 2P = {pi -|- 
P 2 ,Pi G P} is again a lattice polygon. Moreover, in dimension two, any lattice 
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polygon is normal, which means that 2P n = {pi + P 2 ,Pi & P CM?}, which is 
implicitly used in the algorithm. 

Algorithm 3.1. The following algorithm produces a matrix of polynomials in 
(Ti, T 2 , Ta) out of the input polynomials /o,..., /a in variables s = (si, S 2 )- 

• INPUT: A lattice polytope P and polynomials (/o(s)j/i(s),/ 2 (s),/ 3 (s)) 
with no common factor and Newton polytopes Af{fi) contained in P. 

• STEP 1: Consider syzygies {ho,... with JV (hi) C 2 P. Let{h^^\ ... ,h^^^), 
j = \,... ,N, be a 'K-basis of such syzygies. 

(i) (i) 

• STEP 2: Represent the syzygies as linear forms Lj = h]f To+ ■ ■ --{-h f T^. 
Write = X]/ 3 g 2 Pnz 2 ^P/ 3 ®^ switch: 

i P \ i ) 

• OUTPUT: The matrix M of linear forms ■= "^^h^^^Ti. 

We illustrate the steps in Algorithm |3.1| in the following example. 

Example 3.2. Let P be the lattice polygon with vertices (0,0), (0,1), (2,0) and 
(1,1), with lattice points po, •. • ,P 4 , as in the figure. We consider the following four 
polynomials with support in P, where we denote s := (si,S 2 ), and given p := {i,j) 
we write s^ := s^s^.' 

fo = I + 3si + sj + 2 s 2 + S 1 S 2 = sP° + 3sPi + sP^ + 2sP3 + sP\ Ps 
/i = 5sP° - sP^ - sP^ + 2 sP^ - sP\ 

/2 = 7sP» + 3sPi + 2 sP^ + 6sP3 + 3sP% 

/3 = llgPo + OsPi + 4sP" + 3sP3 + 5sP". 

To compute the syzygies in Step 1, we consider the morphism (uq, oi, 02 , aa) 1 —> 
where are polynomials with support in 2P. Let B he the matrix of this 
map in the monomials bases. Since 2P has 12 lattice points and Oifi has support 
in 3P, which has 22 lattice points, then B is a Sylvester matrix of size 22 x 48. It 
can be easily checked that B is full ranked (same as for generic polynomials). Thus, 
the kernel of B has dimension A = 48 — 22 = 26, which is the number of linearly 
independent syzygies. 

To construct the matrix M, assume that we choose as our first syzygy the fol¬ 
lowing A-tuple of polynomials (/ip\ . •., ^a^P with Af {h[^'^) C 2P: 

/jW = _i96s2po + 504sPo+Pi - 257sP«+P2 + 672 sP«+P3 + 234sPo+p% 

/jW ^ _237 sPo+p2 _p420sPo+P3 - 168sP«+p% 

= 28s^P° + 10 sP»+P= - 364sP°+P" + 226 sP'>+p% and 

/jW ^ _2i6sPo+P4_ 

We consider Li = + • • • + h^^pT^ and we write 

Li = (-196To + 28T2)s‘^p° + (504ro)sP«+Pi + (-257To - 237ri + 10 T 2 )sP«+p^ 
+(672ro + 420ri - 364T2 )sP«+P 3 + (234To - I 68 T 1 + 226 T 2 - 216^a)sP“+P^ 
which gives the first column of the 28 x 12-matrix M (computed with Macaulay2 
computer-algebra software |GS| ) 
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-196T_0+28T_2 

504T_0 

-257T_0-237T_1+10T_2 

672T_0+420T_1-364T_2 

234T_0-168T_1+226T_2-216T_3 

0 

0 

0 

0 

0 

0 

0 


0 

-196T_0+28T_2 

504T_0 

0 

672T_0+420T_1-364T_2 

-257T_0-237T_1+10T_2 

234T_0-168T_1+226T_2-216T_3 

0 

0 

0 

0 

0 


0 

0 

-196T_0+28T_2 

0 

0 

504T_0 

672T_0+420T_1-364T_2 

-257T_0-237T_1+10T_2 

234T_0-168T_1+226T_2-216T_ 

0 

0 

0 


12 26 

Matrix (QQ[T , T , T , T ]) <- (QQ[T , T , T , T ]) 

0123 0123 


The columns of this matrix M are given by a choice of a basis of syzygies with 
support in 2P. The corresponding linear forms Li,..., are known as the moving 
planes defining the surface parametrized by fi,..., f^. The associated rational map 
f has deg(/) = 1. It can be checked that the common factor of any maximal minor 
of M is the degree 3 implicit eguation of the closed image of f: 


F = 2643r3 + 2905r2ri + 1345rorf + 91T^ - Sr^Ta - AUT 0 T 1 T 2 + 2 MTIT 2 + • • • , 
as asserted by Theorem \3.3\ below. 

Note that we have to write the lattice points in 2P as a sum of two points in P, 
but in general there is not a unique way of doing this. In our example, for instance, 
Po + P 2 = Pi + Pi, so a choice was made. In fact, it is possible to make a coherent 
choice in general with the use of weight vectors, but any choice will work since in the 
quotient ring A defined in (4.7) below, it holds that XQ.X 2 and X^ are identified. 


We now state our main result. The proof will be given in Section]^ 


Theorem 3.3. Given (/o(s), fi{s), / 2 (s), fsis)) with no common factor, with New¬ 
ton polytopes contained in P and satisfying hypotheses \4.4\ below. Algorithm \3.1\ 
computes a presentation matrix of the implicit equation of the rational map f. 
That is, the rank of the matrix M drops precisely when evaluated at the points in 
the closure of the image of f. 

Moreover, the implicit equation F can be computed as 

(3.1) ^deg(/) _ minors of M). 


The main ingredient for the validity of Algorithm |3.1| to give a matrix repre¬ 
sentation is the choice 2P n of the support of the linear syzygies. Again, the 
“instant” elimination is done in STEP 2, where the s variables give the monomial 
basis which is used to compute the matrix M (and thus they disappear from the 
output!). 

In fact, Algorithm |3.3| can be run without checking Hypotheses |4.4[ We point 
out in Remark |4.7| the possible outcomes. The general algorithm can be refined 
using Theorem 11 in |BDD09] . 


Theorem 3.4. Assume f satisfies the hypotheses of Theorem \ 3.3\ If the lattice 
polygon P can be written as a multiple P = dP' of another another lattice polygon 
P' without interior lattice points, then we can consider in STEP 1 of Algorithm \3. l\ 
syzygies (/iq, ..., h^) with smaller support N{hi) contained in {2d — 1)P' (which is 
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strictly contained in 2P), and the OUTPUT will still be a matrix representation 
for f. Moreover, in case P' is the unit simplex, it is enough to consider syzygies 
with support inside {2d — 2)P'. 


We then have the following comparison between the general syzygy method and 
the naive linear algebra method described in Section 


Remark 3.5. Assume that P is the triangle of size d. Then, as it is enough to 
consider syzygies of degree 2d — 2, they can be found by solving a linear system on 
4(^2^^) variables with {^^) equations. That is, both sizes, as well as the vector space 
dimension of the space of syzygies in this degree, are quadratic in d. The matrix M 
has then a number of rows quadratic in d. The number of its columns equals {'^^^, 
again quadratic in d. Comparing with the sizes in Lemma \2.3\ which are of degree 6 
in d, we observe that the syzygy method is a great improvement on the naive linear 
algebra method! 


The same improvement occurs for any lattice polygon P. Using (2.1), we see that 
syzygies with support in 2P can be obtained by solving a system with approximately 
9vol(P) equations in 16vol(P) variables and both row and column sizes of the 
matrix representation M are of the order of vol(P) and not of its cube, as in 
Lemma 12.31 


3.3. The bihomogeneous case and beyond. As we have mentioned, the main 
motivation for the implicitization problem comes from Computer Aided Geometric 
Design and geometric modeling. In this area, bihomogeneous surfaces (correspond¬ 
ing to rectangular support P) are known as tensor product surfaces, and they play 
a central role, in particular the Bezier surfaces. Quoting Ulrich Dietz [Die98| : “/n 
current CAD systems tensor product surface representations with their rectangular 
structure are a de facto standard”. These surfaces (called NURBS) are given by 
pieces of parametrized surfaces cut by curves. So, it is necessary to control the 
location of the parameter, which can be achieved by computing the kernel of the 
matrix representation we give, as explained in |Busl4j . 

Due to the nature of the base locus of the parametrization, many of the current 
geometric modeling systems do not satisfy the hypotheses to be detailed in |4.4| 
needed for Theorem 3.3 to hold, if considered as homogeneous polynomials (with P 


an equilateral triangle). But if we use a rectangle P as the input in Algorithm |3.1[ 
it is possible to get a full-ranked matrix representation by Theorem |3.3[ In this 
bihomogeneous case, the detailed study of regularity in |BC11] allows to get the 
following improvement in the support of the proposed linear syzygies in STEP 1 
of Algorithm 1 3.1 1 it is enough that the support of these syzygies is contained in a 
polygon obtained by only enlarging the rectangle support P of the input polyno¬ 
mials (approximately) to its double in only the horizontal or the vertical direction, 
instead of considering syzygies with support in (the lattice points of) 2P. 

A more general result can be obtained for bigraded toric surfaces, and in particu¬ 
lar for lattice polygons defining a Hirzebruch surface, that it, for Hirzebruch quadri¬ 
laterals Ha^b.n with vertices (0, 0), (a, 0), (0, b) and {a + nb, b), for any a, bn € N. We 
state this extension in Theorem 13.81 below. 


Algorithm 3.6. Take as INPUT a Hirzebruch lattice polygon P = Ha^b,n o,nd 
bivariate polynomials (/o, /i, /2, /a) with Newton polygons contained in Ha^b,n> o-nd 
which satisfy the hypotheses\4.8[ Run algorithm\3. 1\ with the following modification: 
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in STEP 1 consider a basis of syzygies with support in the smaller lattiee quadri¬ 
laterals H 2 a-i,b-i.n (or Ha-i^ 2 b-i,n instead). T/ie OUTPUT is the corresponding 
matrix M of linear forms. 

In most cases, it is convenient to consider syzygies with support in H 2 a-i,b-i,n 
rather than in since the first one has less lattice points. 

{ny-\-x,y) 

In general, a Hirzebruch lattice polygon Hx^y^n ’ 
has the shape in the diagram on the right. 

( 0 , 0 ) (®, 0 ) 

Remark 3.7. Note that for n = 0, Ra,6,o is the standard lattiee rectangle with 
vertiees in (0, 0), (a, 0), (0, b), (a, b), and thus Algorithm 3.6 works in particular in 
a standard bihomogeneous setting. 

Theorem 3.8. Given (/o(s), /i(s), f 2 {s), fd,{s)) with no common factor, with New¬ 



ton polytopes contained in a Ha^b,n ond satisfying hypotheses 4-8 below, Algorithm 3.6 
computes a presentation matrix of the implicit equation of the rational map f. That 
is, the rank of the output matrix M drops precisely when evaluated at the points in 
the closure of the image of f. 

Moreover, the implicit equation F can be computed as 

(3.2) ^deg(/) _ gQ(j^)jji(iximal minors of M). 

The proof of Theorem |3.8| will be also given in section 


4. The hypotheses via toric geometry and the proofs of our main 

RESULTS 

In this section we will recall a minimum of theoretical tools from toric geometry 


in order to state the hypotheses needed for Theorems |3.3| and |3.8| and to give their 
proofs. The main homological commutative algebra tools that are the core of the 
proofs are recalled in Appendix [X| 

We refer to |Cox03bllFul93[[CLSlI| and [GKZ941 Ch.5&6] for the general notions, 
and to |KD061 §2], |BDD09llBotlla| for applications to the implicitization problem. 
Any reader only interested in the application of Algorithm |3 .1 1 or its bihomogeneous 
(toric) refinement given in Algorithm |3.6| can skip this section. 

As usual, we denote by K* = IK\ {0} the multiplicative group of units of K. The 


first observation is that we can equivalently consider our parametrization ( 1 . 1 ) as 

(IK*)2 - 


a map / : —+ P'^(K) or / : (K*)^ P'^(K) with image inside 3-dimensional 

projective space, and domain a dense open set U in affine space or the torus 
(K*)^ over K, given by 


(4.1) 


S ^ (/o(s) : /i(s) : / 2 (s) : fsis)), 


for any s € U, and we have the commutative diagram 
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In fact, if F is the implicit equation of the (closure of the) image of /, the (closure 
of the) image of / is the closure of under the standard embedding P^(K) 

and its equation is the homogenization of F. 

Similarly, we can consider our rational parametrization from any algebraic vari¬ 
ety which contains the domain of / as a dense subset. We will choose embedded 
or abstract compact toric varieties to get a degree or multidegree notion that will 
allow us to get homological arguments to “bound” the support of the syzygies in 
Theorem 18.31 and in 18.81 


4.1. Toric embeddings. Let / be a rational map as in (4.1). The base points of 
the parametrization are the common zeros of /o,..., /a, that is, the points where 
the map is not defined. We assume that / is a generically finite map onto its 
image and hence it parametrizes an irreducible surface ^ C P^. We also assume 
without loss of generality that gcd(/o,/i,/ 2 ,/a) = 1, which means that there are 
only finitely many base points. 

Let P C be a lattice polygon with m -|- 1 lattice points, which contains the 
Newton polygons Af(fo), ■ ■ ■, A/'(/a). Write PnZ^ = {po,... ,Pm}. The polygon P 
determines a projective toric surface C P™ as the closed image of the embedding 


(IK*)2 

(si,S2) 


l-A- 


...) 


where i = 0,... ,m. For example, the unit triangle with vertices (0,1), (1,0) and 
(0,0) (or any lattice translate of it) corresponds to P^, and any lattice rectangle 
gives a Segre-Veronese projective embedding of P^ x P^, which are special cases of 
toric embeddings. 


Example 4.1. Assume P is the unit square, with m -I- 1 = 4 integer points: 

Po = ( 0 , 0 ), Pi = (1,0),P2 = (0,1),P3 = (1,1)- 

A polynomial fi with Newton polytope contained in P looks like 
(4.3) fi{s) = 0(0,0) + 0(1,0)51 + 0(0,1)52 + a(l,l)SiS2. 

We take 4 new variables (Vq : Xi : X 2 : X 3 ) as the homogeneous coordinates 
in P^. The toric variety is the projective variety in P^ cut out by the relation 
XqXo — X 1 X 2 = 0 . This binomial equation comes from the primitive affine relation 
P0+P3 = Pi+P2, which implies the multiplicative relation sP°sP^ = s^^s^^ between 
the monomials with these exponents. The coordinate ring of !Xp is the quotient ring 
K[Vo,...,V3]/(XoX3-V 1 V 2 ). 

In general, we will call (Vq : • • • : X^) the homogeneous coordinate variables in 
P™. Write P n = {po, ■ • ■ ,Pm}- We set one variable Xi for each integer point 
Pi in P and we record multiplicatively (by binomial equations) the affine relations 
among these points. These binomials generate the toric ideal Jp = J(f7p), which 
defines the variety 3i^p C P™. To each 

m 

(4-4) fi{s) ='^apy', 

we associate the homogeneous linear form 

m 

(4.5) gi{s) ^y^^ap.Xj. 

2=0 
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For instance, in Example |4.1[ the polynomial fi gets translated to 

gi{Xo, ..., X 3 ) = a(o,o)-^o + a(i,o)-’^i + <^{o,i)X 2 + a(i,i)-^3) 
and over we have the relation X 0 X 3 — X 1 X 2 = 0. 

The rational map / factorizes through .^p in the following way 

(4.6) (]k*) 2^p3 

P < 

y 

P / 

where g = {go ■ gi ■ g 2 ■ 93 ) given by four homogeneous linear polynomials 
go, gi, 92,93 in m+ 1 variables. Thus, we have a new homogeneous parametrization 
g of the closed image of / from ,!^p. The polynomials gi generate an ideal I in the 
coordinate ring 

(4.7) A = K[Xo,...,Xrn]/Jp 

of ^p. This ideal / defines the structure of the base locus in ^p. 

The embedding p : (K*)^ —>■ provides a Z-grading in the coordinate ring A of 

3^P, which is used to study the map g with the tools recalled in Appendix [X} 

4.2. Abstract toric varieties and Cox rings. Given a lattice polygon P, one 
can also associate to it an abstract compact toric variety Xp that naturally contains 
the torus (K*)^ as a dense open set (via the map we call j below), adjoining a torus 
invariant divisor to each edge of P. We refer the reader to |CLS111 ICL098| for the 
theory and details. 

The map / also defines a rational map / that makes the following diagram 
commutative: 

(4.8) Xp 

A \ — 

'' N / 

J \ 

U A 

(K*)2 - - ^P3 

P / ^ 

/ 

P / 

^P 


Example 4 . 2 . Assume P is the unit square, with N = A edges: the segments Ei = 
[(0,0), (1, 0)], E2 = [(0,0), (0,1)], Eo = [(0,1), (1,1)], and E^ = [(1, 0), (1,1)]. The 
respective inner normal vectors r]i = (0,1), 772 = (1,0), 773 = (0,-1), 774 = (—1,0) 
satisfy the linear relations 771 + 773 = 0, 772 + 774 = 0, which give rise to two homo¬ 
geneities. We introduce four associated variables Y = (Yi, ..., F4). A polynomial fi 
with Newton polytope P as in ( 4 . 3 ) defines a hihomogeneous polynomial (in (Yi, Y3) 
and (^2,^4 )): 

fi{Y) = a(o,o)^^4 + <2(1, 0)^1 ^4 + a(o,i)Y2Y3 + a(i_i)Yiy2- 
These polynomials fi define the map f = (/g : ..., /g). 


The main motivation for this change of perspective comes again from the com¬ 
mutative algebra results needed for the proof of Theorems |3.3| and |3.6[ The Cox 
ring of Ap is endowed with a more natural multigrading, which is finer than the 
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grading obtained via the embedded projective variety Also, this point of view 
has an impact in the computations, as the number of variables to eliminate is 
smaller (one for each edge of P, instead of one for each lattice point in P). In our 
small example |4.2[ there are four edges and four lattice points, but the number of 
edges can remain constant while the number of lattice points goes to infinity. 


4.3. Precise hypotheses and proof of Theorem | 3 . 3 [ In this subsection we 
detail the precise hypotheses that ensure the validity of Theorem |3.3| and we prove 
it, based on results in |BDD09| . We first need to recall a few standard definitions 
from commutative algebra. 


Definition 4 . 3 . Given (nonzero) homogeneous polynomials {go, ■ ■ ■ ,gs,), defining 
a rational map g : I 

is a common zero set of go, 


as in (4.61, a point p G iTp is a base point of g if it 
,go, that is, if p is a zero of the ideal I C A in 3"p. 


Let p G be a base point of g. The local ring of p is the ring Ap = {hi/h 2 , hi G 
A, /i 2 (p) ^ 0}, with the natural operations induced from A (in turn, naturally in¬ 
duced from the polynomial ring). Let Ip be the ideal generated by (the classes of) 
go, ■ ■ ■ go in Ap. We say that p is a loeal complete intersection base point if Ip can 
be generated by only 2 elements. We say that p is an almost complete intersection 
base point if Ip can be generated with 3 elements. 

We have similar definitions for the map f : (K*)^ ---» P^. 


For a given lattice polygon P, here are the hypotheses we need in terms of g: 


Hypotheses 4 . 4 . There are only finitely many base points of g on 3Lp which are 
local complete intersections. 


We cannot easily find hypotheses on / equivalent to Hypotheses |4.4[ Given a 
lattice polygon P, an edge E of P, and a polynomial fi with N{fi) contained in P 
as in (4.41, the restriction of fi to E is defined as the sub-sum of the monomials 
with exponents pi in E. We have the following partial translation. 


Proposition 4 . 5 . Let f,l7' and g be as in (4.6|. Then 


(1) There are only finitely many base points of g on SLp if and only if there are 
only finitely many isolated base points of f in the torus and for each edge 
E of P, at least one of the restrietions is nonzero. 

(2) If g has finitely many isolated base points on STp whieh are local complete 
intersections, then the base points of f in the torus are local complete in¬ 
tersections. 


Proof. The map p defines an isomorphism between (K*)^ and its image (which is 
an open dense subset of i^p), sending a base point q of f (that is, a point where 
foiq) = ‘ ‘ ‘ = fsiq) = 0) to a base point p = p{q) of g, and reciprocally, any base 
point p of g in p(IK)^ is the image of a base point q of /. Moreover, we have an 
isomorphism between the ideal generated by /o,..., /a at q and Ip. Any base point 
of g outside the image of p cannot be seen in the torus. But these points are either 
the fixed torus points corresponding to the finitely many vertices of P, or they lie 
in the torus of the toric divisor Dp fn. STp associated to an edge E oi P. As Dp has 
dimension 1, there are finitely many solutions as long as at least one of the fi^p is 
nonzero. □ 


The following example shows that the converse 
does not hold. 


of item in Proposition |4.5 
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Example 4.6. Consider the parametrization with 6 monomials: (fo, fi, f 2 , fs) = 
{sC + 2, sC — 3st^, st"^ + 5s'^t^,2 + s'^t^). Then, f has no base points in the torus. But 
if we consider their standard homogenizations to degree 8 polynomials (that is, we 
take P equal to 8 times the standard unit simplex in the plane), the corresponding 
homogeneous polynomials g^,... ,g^ have one base point “at infinity” which is not 
even an almost locally complete intersection. 

We now give the proof of Theorem |3.3| 

Proof of Theorem \3.3[ Given (/o(s), fi{s), f 2 {s), fsis)) with no common factor and 
Newton polytopes contained in P, the corresponding polynomials gi associated to 
fi are homogeneous of degree d = 1 and satisfy Hypotheses |4.4[ 

From |BDD09[ Cor. 14] one has that for d = 1, the matrix of coefficients of a 
if-basis of the module of Syzygies of g in any degree v >2 with respect to a iF-basis 
of the graded piece A^, of A, is a matrix representation for the closure of the image 
of /, which equals the closure of the image of g. 

In particular, we can take v = 2. In STEP 1 of Algorithm |3.1[ the syzygies 
{hQ \ ..., h^'^) ioT j = 1,..., N with A/’(/ip^) C 2P for all i,j, provide a K-basis of 
the module of syzygies of g in degree 2, since classes of monomials of degree 2 in A 
correspond to monomials in the s variables with exponents in 2P. 

Equality ( |3.1[ ) follows from Theorem 13 in [BDD09] . □ 

In principle, given a rational map /, we could take any lattice polygon P con¬ 
taining the union N'{f) of Newton polytopes of /o, -.., /a- Note that the hypothesis 
that / is generically finite implies that Af{f) is two-dimensional. Taking P strictly 
containing Af{f) will increase the number of exponents and will in general produce 
bad behaviour of g at the fixed points in corresponding to the vertices of P 
which do not lie in JV{f). 

Remark 4.7. We can check algorithmically if fo,..., fs, have finitely many so¬ 
lutions over (K fi and if for any edge E at least one of the restrictions is 
nonzero. So, by Proposition \4.^ we can check whether g has finitely many base 
points in A^p. 

Assume the dimension of the base locus of g is zero. As we remarked in Exam- 
pie\r^ even if we could check the local behavior of the base points of f in the torus, 
this would not imply the satisfiability of Hypotheses But what if we don’t check 
this and run Algorithm \3.1\ ? ... nothing bad! 

We then check whether the output matrix M in has full rank. 

• If the rank of M is not maximal, then there is at least one base point p of g 
which is not an almost local complete intersection. In this case, we cannot 
get the implicit equation, but we get a certificate of the bad behavior of the 
base locus (without computing it). 

• If the rank of M is maximal, it may happen that the its rank drops when 
evaluated at points outside due to the existence of an almost complete in¬ 
tersection but non complete intersection base point. In this case, the greatest 
common divisor of the maximal minors of M would have irreducible factors 
other than the implicit equation E. In fact, the existence of other irreducible 
factors is equivalent to the fact that there exists a base point which is an 
almost local complete intersection but not a local complete intersection. 
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4.4. The hypotheses and proof of Theorem |3.8[ In this subsection we detail 
the precise hypotheses that ensure the validity of Theorem |3.8| and we prove it, 
based on results in [Botllbj . 

Given P and /, here are the hypotheses we need in terms of the map / in 


(4.8): 


Hypotheses 4.8. There are only finitely many base points of f on Xp which are 
local complete intersections. 


Again, we cannot easily find hypotheses on / equivalent to Hypotheses |4.4[ since 
good algebraic behaviour of the base points in the torus does not imply the same 
behaviour for the possible base points of / at the invariant divisors in Xp associated 
with the edges of P. 


Proposition 4.9. Let f,Xp and / as in 4-8 Then 

(1) There are only finitely many isolated base points of f on Xp if and only if 
there are only finitely many isolated base points of f in the torus and for 
each edge of P, at least one of the restrictions of the fi is nonzero. 

(2) If f has finitely many base points on Xp which are local complete intersec¬ 
tions, then the base points of f in the torus are local complete intersections. 


We next give the proof of Theorem |3.8 


Proof of Theorem \ 3 .^ By hypothesis, there are only finitely many isolated base 
points of / on the toric variety Xp associated with P := Ha^b,m which are local 
complete intersections. There are four primitive inner normal vectors of P: 771 = 
(0,1),772 = (0,1),773 = (—1,0),774 = (—l,n), which satisfy the linear relations 
773 = —771,774 = 77772 — 771. So any multidegree v can be described by a “bidegree” 
(1^11 ^2) given by the degrees with respect to the first normals and which fixes (up 
to translation) the associated polytope Pi, with the same normals as P. Thus, by 
[Botllbl Thm. 5 . 5 ] the matrix of coefficients of a AT-basis of the module of Syzygies 
of / in any bidegree {vi, V2) with vi > 2a — 1 and V2 b — 1 (or vi > a — \,i'2 > 
2b — 1 ). with respect to a X-basis of the bigraded piece (z/i, 7^2) of the Cox ring 
of Xp, is a matrix representation for the closure of the image IX oi f (which equals 
the closure of the image of /). 

Taking (pi,P 2 ) = (2a — 1,6 — 1) one has that in STEP 1 a basis of syzygies 
(6 ,q^\ ..., for j = 1,..., X with all J\f{h[^'^) with support in the quadrilateral 
H 2 a-i,b-i,n, provides a K-basis of the module of syzygies of / in bidegree (2a — 
1,6— 1). Hence, the matrix M of coefficients of such syzygies obtained in STEP 2 
gives a representation matrix for 

□ 


Equality (3.2) also follows from Theorem 5.5 in [Botllbj . 


5. Examples 


This section consists of four examples which highlight the usefulness of our ap¬ 
proach. Example 5.1 is taken from a case studied in [BDD091 Ex. 18] of a sparse 
parametrization where projective implicitization does not work due to the nature 
of the base locus of the map, but Algorithm |3.1| is applicable with a right choice of 
polygon P read from the monomials of the input polynomials. In Example |5.2| we 


^The choice of the bidegree is less obvious than in the graded case. For further details, see 
definition of Dtp ( 7 ) in IBotllbl Thm. 5.5], or the analysis of the bidegree in the standard bigraded 
case in [Botllbl Sec. 7.1] 
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show how the method in Algorithm |3.1| works for a parametrization given by fewno- 
mials of high degree, where classical resultant tools fail due to the computational 
complexity. In Example |5.3[ classical resultant tools fail because of the existence 


of a base point in the torus. Finally, in Example 5.4 we compare the methods in 


Algorithms |3.1| and 3.6 


5.1. A very sparse parametrization. Consider the parametrization with 6 mono¬ 
mials given in |BDD09l Ex. 18]: (/o,/i,/ 2 ,/a) = + 2, st^ — 3st^, st'^ + ,2 + 

s^t^). The matrix representation can be computed using the package MatrixRep- 
Toric.m2 |BD10| in the computer algebra software Macaulay2 |CS| . 

One first defines the map / given by polynomials in the ring S = Q[s,t] (note 
that for easiness of typing, we call the variables (s,t) instead of (si,S 2 )): 

S = qq[s,t] ; 

f = {s*t"6+2, s*t~5-3*s*t~3, s*t"4+5*s"2*t"6, 2+s~2*t~6}; 

Consider P the lattice triangle with vertices 
(0,0), (1,6) and (2,6). 

One can compute P by the command: 

P = polynomialsToPolytope L 

The lattice-points of P can be computed using the 
auxiliary Macaulay2 package Polyhedra as: 
latticePoints P 

By taking syzygies with support in 2P, on gets a matrix representation of size 
17 X 34. The greatest common divisor of the 17-minors of this matrix is the homo¬ 
geneous implicit equation of the surface: 

2809ro^T(‘ -P 124002Tf - 5618ro^Tf T 2 -P 66816ToT]^r2 -P 2809TffT| 
-50580ro^Tf r| -P 86976T(‘r| -P 212T^T^ - 14210ToTf -P 
+ imi2TlT^ + lieToTl -P 841T| -P IdOdSTo^T^Tg - 169849Tor]^r3 
-14045ro‘T2T3 -P 2Q1^21T^TIT2T3 - 468288T^T2T3 - 12Q8T^TITj, 
-Pl57155ToTf r|r 3 - 8 im 8 T^TlT 3 - 129215Tf T|r 3 - 4528ToT2^r3 
-12673r|T3 - 16695ro^Tfr| -P 169600T]^r| -P 30740ro^T2T| 
-A3,888^TqTIT2TI -P 82A3,AT^tIt^ + 269745T^TiT^ + 36696ToT^T^ 
-P63946T2‘7’3 + 2775ToT^T^ - 19470To2T2r3^ -P 177675r2T2T| 
-85360roT|?1* - 109490T|r| - 125T^T^ + 2900ToT2T^ 

+7325TiT^ - I 25 T 2 TI 

We can set Tq = 1 to get the affine equation. 

The map g is computed with the following command: 

g = teToricRationalMap f; 

The matrix representation and the implicit equation are computed as follows: 

M = representationMatrix (teToricRationalMap f,2); 
implicitEq (L,2) 
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Notice that the 2 as second parameter in the computation of M is precisely the 2 
in the support 2P of the syzygies. For a deeper understanding of the choice of the 
parameter 2, see Appendix [X} 

In the language of Section the coordinate ring associated to is A = 
K[Xo,...,X 5 ]/Jp, where Jp = (X| - X 2 Xi,X 2 X 3 - X^X^^X^ - XiX 3 ,Xf - 
XqXz). The parametrization g over ^p is given by {go, gi, g 2 , 93 ) = {^Xq + 
X^,2Xq + X 4 ,— 3 X 1 + X^,X 2 + 5 X 5 ). This matrix can be computed as the 
right-most map of the pq = 2 d = 2 strand of a graded complex as explained 
in Appendix The method fails over (i.e P = the triangle with vertices 
(0, 0), ( 8 ,0), (0, 8 )) due to the nature of the base locus. One can see this just by 
computing a matrix representation and verifying that it is not full-ranked. 

5.2. Fewnomials with high degree. This example contributes to show how the 
method works fine for high degree fewnomials involved in the parametrization. 

Consider the polynomials (/o, /i, / 2 , /a) = (l-l-st-|-s^'^, s'^+s^'^, s®^). Let 

/ : be the parametrization that maps {s,t) >-)• (/ 1 // 0 ,/ 2 // 0 ,/ 3 //o)(s, t)- 

The implicit equation of the closure of the image of / could be computed by elimi¬ 
nating the variables (s, t) as follows (using general elimination procedures based on 
Grobner bases in Macaulay2): 

R = QQ[s,t, X, y, z, w] 

fO = 1 + s*t + s~37; fl = s~7 + s~47; f2 = s~37 + s~59; f3 = s~61 
LI = x*fl - y*f0; L2 = x*f2 - z*f0; L3 = x*f3 - w*fO 
eliminate ({s,t}, ideal(Ll,L2,L3)) 

In a 2014 standard desktop computer this routine does not end after one hour of 
computation. We also tried the well implemented eliminate command in Singular 
[DGPSI 2 ] , but with the same lack of answer after a couple of hours of computation. 

By homogenizing with an auxiliary variable u we could try eliminate {s,t,u) 
using Macaulay resultant methods, but we easily figure out that the homogenized 
forms Li,L 2 ,L^ vanish identically over the point {s,t,u) = (0:1:0). This implies 
in particular that the Macaulay resultant IleS(^s,t,u){Li, L 2 , L 3 ) is identically zero. 

Finally one can compute the implicit equation (and matrix representation) by 
implementing the implicitization techniques described in this article. A not very 
efficient (but efficient enough) routine in |Botl0| gives the toric map g in less than 
2 minutes and the desired matrix representation M in less than one more minute. 


5.3. Fewnomials with base points in the torus. This example shows a case 
where classical resultants cannot be applied to compute the implicit equation, but 
the techniques in the paper can. Anyway, we recall that the aim of the matrix 
representations is to provide a better and more complete tool for representing a 
surface, and hence, the point presented with this example is just one extra advantage 
of the method. 

Consider the following parametrization (/o, /i, / 2 , /a) = (1—is, —ts^^+1, —i(—s^®-|- 
— t) given by four polynomials that define a parametrization / : 
that maps (s,i) 1 —)■ (/ 1 // 0 ,/ 2 // 0 ,/ 3 //o)(s, i). The implicit equation cannot be 
computed by eliminating the variables (s, t) with classical resultants, because the 
point (1,1) is in the base locus. However, this fact does not imply any problem in 
Algorithm |3.1| 


With Algorithm 3.6 with a rectangular P = 7 / 38 , 2 , 0 , it takes 0.058 seconds in a 
standard 2014 desktop computer to get a matrix representation M of the closure 
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of the image of /. The size of M is 152 x 194 and the gcd of its maximal minors 
has degree 110 . 


5.4. Comparison with and without embedding. While Algorithm |3.1| holds 
with great generality, when dealing with polynomials with rectangular support 


(which can be interpreted as bihomogeneous polynomials), Algorithm 3.6 provides 
a smaller matrix representation. 

Consider the following four polynomials fo, ■ ■ ■, fs- 
fo = 3sf S 2 - 2siS 2 - S? + SlS 2 - 3si - S 2 + 4 - S 2 


/l — 3 s^S2 — sf — 3siS 2 — Si + S2 + ^2 + S 2 4” 


*1*2j 


/2 = 2 sjS 2 — 3s'fs2 — s'f + S 1 S 2 + 3si — 3 s2 + 2 — S 2 , and 
/a = 2sfs| — 3sfs2 — 2 siS 2 + sf + 5siS2 — 3si — 3s2 + 4 — s^. 

The Newton polytope F = A/’(/) is the rectangle {(x,y) : 0 < x,y < 2)}. If 
we apply Algorithm 3.1 (as we illustrated in Example 3.2), we obtain a matrix 
representation M of size 25 x 51. 

The associated toric variety Xp can be identified with the (2, 2) Segre-Veronese 
embedding of x P^ in P ^ (se e |BD07l IBDDOQI iBotllaj l. 

By means of Algorithm |3.6| we get a matrix representation M from a basis of 
linear syzygies of bidegree (2.2 — 1,2—1) = (3,1). This matrix representation can 
be computed using the algorithm developed in [BotlO] and implemented in M2, as 
the matrix M^, for bidegree v = (3,1), and one obtains a square 8 x 8-matrix. Its 
determinant equals the implicit equation F\ 


F = 63569053X - 

0 


159051916X X 
0 1 


6 2 

175350068X X 
0 1 


5 3 

82733240X X 
0 1 


4 4 

2363584X X 
0 1 


Notice that the matrix Al (3 1 ) is considerably smaller than the 25 x 51-matrix M 
because instead of considering syzygies with support in the rectangle 2P = {(a;, y) : 
0 < a;, ?/ < 4}, the syzygies are taken with support in the smaller rectangle {(a;, y) : 
0 < a; < 3, 0 < 2 / < 1}. 


A. Appendix: Commutative algebra tools 

This appendix is devoted to highlight the tools of homological commutative alge¬ 
bra and algebraic geometry that are needed to justify the validity of Algorithms |3.1| 
and |3.6[ and to explain the choice of the support of the syzygies in STEP 2 which 
define the matrix representation of the parametrized surface. 

Given P, the toric embedding p : (K*)^ —>■ P"* in Section 4.1 provides a Z- 
grading in the coordinate ring A = IK[Ao,..., Am]/Jp in (4.7) of which can 
be used to study the map g in Diagram |4.6| and its associated Rees and symmetric 
algebras, denoted by Rees^(/) and Sym^(/) respectively. Notice also that the 
graded ring A coincides with the affine semigroup ring of the lattice polytope P, 
which is Cohen-Macaulay and normal because P has dimension 2. 

The grading in A plays a key role in the elimination process. The matrix rep¬ 
resentation M of Section |3.2| depends on a choice of degree u, as was shown in 
Section |3.1| for the case of curves. The reason why u needs to be considered is 
rather technical, and a complete explanation involves sheaf cohomology. From a 
more naive point of view, the implicit equation of the surface := ixa^g) is written 
in the variables T = (Tq, ..., T 3 ) but depends on the algebraic relations among the 
polynomials gi, which lie in A. Fixing a degree v in A can be thought as eliminating 
the variables of A, by hiding them in the monomial basis of the graded piece Ay. 
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In turn, recall that the variables X = (Xq, ... ,Xm) in A code monomials in the 
original s variables, with exponents in the lattice points in P. 

More geometrically, consider the graph variety F of 5 where both group of vari¬ 
ables X and T are involved. The elimination process can be understood geometri¬ 
cally as projecting F via 1^2 in the following diagram 


F C ^ X 

TTl 

3^p -= ini(/)- 

In the correspondence between subvarieties of xP^ and bigraded algebras, the 
inclusion of the graph F C x P^ corresponds to the surjection A\Tq, Ti, T 2 , Ts] ^ 
Rees^(/), the Rees algebra of the ideal I generated by ... ,g 3 over the coordi¬ 
nate ring A. The projection 712 (F) corresponds to eliminating the variables Xi of 
ReesA(F). We denote by 1(712 (F)) the defining ideal of 712 (F) C P^. 

How to eliminate the X variables from ReesA(F) algebraically? A standard pro¬ 
cedure is to find a free graded presentation Fi F 2 —> Rees^(/) —>■ 0 and a 

degree p (in the X variables) such that the Fitting ideal generated by the 

maximal minors of M^, (in the graded strand {Fi)^ {^ 2 ) 1 ' —t ReesA(/)i, —>■ 0 ) 

computes 1(712 (F)). It happens that no universal way to compute such a free pre¬ 
sentation is available, so the idea is to “approximate” Rees^(/) by the (hopefully) 
similar algebra Sym^(/) that admits such a universal resolution. These resolutions 
of the symmetric algebras are known as approximation complexes, they were intro¬ 
duced in [HSV821 IHSV83a] and their application on elimination theory was done 
in |BJ031 IBus06| . The last map of the approximation complex is the following in 
our case: 

Zi[ro,ri,T2,r3] ^ A[ro,ri,T2,r3] ^ Sym^(j) ^ o, 

where = {( 00 , 01 , 02 , 03 ) S A^ : 'Y^aigi = 0 } is the first module of syzygies of 
90 , 91 , 92,93 and M' = [Tq Ti T 2 T 3 Y, that is, 

M' ■ (oo, 01 , 02 , 03 ) := ^ ttiTi. 



The cokernel of M' is Sym^(/) = A[To, Ti, T 2 , T 3 ]/J, where 

J := {^OiTi : Oj G A[To,Ti,r 2 ,T 3 ] and '^atgi = 0}. 

We can recognize the origin of the linear forms Li in STEP 2 of our algorithms! 

But there is a remaining question: which is the relation between Reesyi(/) and 
Sym^(/)? Which variety does Sym^(J) define? Can we use 5^(M') to compute 
^ 71-2 (r) for some p7 The answer is that in case there are finitely many base points 
and for each base point p, the local Ip is a local complete intersection, then ReesA(F) 
and Sym^(/) define the same scheme in .^x (thus, in xP^). As Rees^(/) is m- 
torsion free, both algebras coincide module the m-torsion of Sym^(/), Rees a{I) — 
Sym^(/)/iF^(Sym^(/)), where m is the maximal ideal generated by Xq, ... ,Xm- 
Thus if / is a local complete intersection and ly is such that 7F^(Sym^(/));^ = 0, then 
ReesA{I)u — Sym^(/)y. This happens for p > vq := 2 hy Theorem 11 in [BDD09] 
(in fact, that result also proves Theorem 3.4 as remarked before). In particular, in 
this case, (?(M(,) computes T 7 r 2 (r) for any p > vq. 
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In fact, the condition of I being a local complete intersection can be relaxed to the 
condition of being locally an almost complete intersection, (i.e. Ip can be generated 
by 3 elements, for any p in the finite set V{I)). In this case, dim(Sym^(/)) = 
dim(Rees. 4 (/)). Since there is always a surjective map Sym^(/) ^ Rees^(/) then 
R(Sym^(/)) = r U [/, where U has the same dimension. In particular, 7r2(r) U 
7 r2(C/)7r2(y(Sym^(/)). For i/ > i/q, Sym^(J) is m-torsion free, and S'(M') computes 
/(7r2(R(Sym^(J))) for any ly > vq. So, the gcd H of the maximal minors of MV 
contains has the homogenization of the implicit equation V as a factor. 

In the bigraded case of Hirzebruch surfaces, in particular in the standard bi- 
graded case, the basic ideas are similar, but new technical details have to be man¬ 
aged in order to determine the bidegrees for which the torsion of the symmetric 
algebra vanishes. We refer the reader to |Botllb| for the details and proofs. 
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